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Abstract 

We  show  that,  owing  to  small  noise  and  smoothed  excursions,  the  unbiased  FIR  filtering 
algorithm  produces  more  accurate  estimates  of  the  TIE  model  of  a  local  clock  than  the 
sawtooth-less  GPS-based  measurements.  The  algorithm  becomes  much  more  efficient  in 
estimating  the  fractional  frequency  offset  and  higher  order  states.  For  the  algorithm  to 
work  with  minimum  errors,  the  optimal  time  step  Topt  and  number  of  the  points  in  the 
average  Nopt  must  be  ascertained  in  the  sense  of  the  minimum  mean-square  error  (MSE)  for 
each  of  the  clock  states.  The  relevant  dependencies  of  Nopt  on  Topt  are  found  experimentally 
for  a  crystal  clock.  The  trade-off  between  the  estimates  and  sawtooth-less  measurements  is 
illustrated  graphically. 


INTRODUCTION 

Evaluation  of  the  time  interval  error  (TIE)  of  a  local  clock  is  easily  provided  using  the  commercially 
available  GPS  timing  receivers.  The  problem  is  that  the  GPS-based  measurements  with,  for  example,  the 
Motorola  family  of  the  receivers  such  as  the  SynPaQ  III  GPS  Sensor,  induce  the  sawtooth  noise  of  about 
±50  ns  owing  to  the  principle  of  the  one  pulse  per  second  (1PPS)  signal  formation  utilized  in  the 
receiver.1  By  means  of  the  receiver  [1]  or  auxiliary  software  [2],  the  sawtooth  noise  may  be  reduced  by  a 
factor  of  about  10.  In  the  measurements  of  the  TIE,  it  can  also  be  suppressed  using  the  recently  designed 
unbiased  finite  impulse  response  (FIR)  filter  and  algorithm  [3,4]  allowing  for  an  error  less  than  that 
obtained  by  a  standard  Kalman  filter.  The  first  comparative  results  in  estimating  the  TIE  have  shown 
[5,6]  that  the  FIR  estimates  become  much  more  efficient  than  the  sawtooth  correction  if  the  unbiased  FIR 
filter  is  optimized  for  the  time  step  Topt  (obtained  by  thinning  in  time  the  1PPS  pulses)  and  the  number  of 
the  points  Nopt  in  the  average. 

In  this  paper,  we  take  a  more  precise  look  at  the  problem  comparing  the  TIE  estimates  of  a  local  crystal 
clock  to  provided  simultaneously  the  GPS-based  sawtooth-less  measurements  and  measurements  for  a 
reference  rubidium  clock.  We  also  present  the  results  of  the  experimental  evaluation  of  the  optimal  Nopt 
and  Topt  for  several  states  of  the  local  clock,  thereby  answering  the  question:  what  are  optimum  values  iopt 
and  Nopt  for  the  unbiased  FIR  filtering  algorithm  in  the  sense  of  the  minimum  mean-square  error  (MSE)  in 
each  of  the  estimates? 


!Clark  and  Hambly  reported  in  [1]  that,  in  the  modified  receivers,  the  noise  bounds  are  reduced  to  ±13  ns. 
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AN  UNBIASED  FIR  FILTERING  ALGORITHM 

Before  discussing  the  results  of  the  experimental  studies,  it  is  in  order  first  to  describe  in  brief  the  TIE 
model  of  a  local  clock  and  the  unbiased  FIR  filtering  algorithm  proposed  in  [4]  and  thereafter  studied  in 

[6]. 

TIE  Model  of  a  Clock 

Most  commonly,  the  TIE  polynomial  model  of  a  clock  projects  ahead  on  a  horizon  of  N  points  from  the 
starting  point  n  =  0  with  the  Kth-degree  Taylor  polynomial: 

k  TPnP 

xx  (n)  =  I  x  , - —  +  w,  (n,  T) 

P= 0  p\ 

-  xx  +  x2t  n  +  —  r2n2  +  —  T2n2  ...  +  wx{n,z),  (1) 

2  6 

where  X/+  \  =  x/+  \  (0),  /  e  [0,  K\,  are  initial  states  of  the  clock  and  w\{n,x)  is  a  clock  noise  with  known 
properties.  By  extending  the  time  derivatives  of  the  TIE  model  to  the  Taylor  series,  the  signal  and 
observation  equations  become,  respectively, 

A(n)  =A(n)A(0)  +  w  (n,z),  (2) 

g(n)  =  CA(n )  +  v(«),  (3) 

where 

X(n)  =  | \xx{n)x2(n)  ...  xK+l(n)]T 
is  a  (K  +  1)  x  l  vector  of  the  clock  states  and 

1  Tn  T2n2  / 2  ■■■  ( Tn)K/K\ 

0  1  Tn  •••  (rn)K-1 

A(n)  =  0  0  1  •••  (Tnf-2 /(K- 2)1  .  (4) 

0  0  0  •••  1 

is  a  (K  +  1)  x  (K  +  1)  time-varying  system  matrix. 

For  an  equal  number  of  the  states  and  measurements,  the  observation  vector  is  described  by 

£(»)  =  [<?,(»)  &  (»)-  4-,  (n,r)f 

and  the  (AT+l)x(AT+l)  measurement  matrix  C  is  typically  unit.  The  {K  +  1)  x  1  clock  noise  vector  is 
described  by 
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w  (n,z)  =  [wx(n9T)  w2(n,T)...  wK+l(n,r)]T 
with  the  components  caused  by  the  oscillator  noises.  Finally,  the  noise  vector 

\{n)  =  [vl(n)v2  («)...  vK_x{n)\T 


contains  correlated  or  uncorrelated  components  that  are  not  obligatory  white  Gaussian.  The  GPS  noise 
v  (n)  dominates  on  a  horizon  N ;  that  is,  typically,  {w2u  (n,T))N  «  ( v ]  ( n))N .  Therefore,  w  ( n ,  x)  is 
neglected  in  (3)  and  in  the  FIR  procedure  [4]. 

An  Unbiased  FIR  Filtering  Algorithm 


The  unbiased  FIR  filtering  algorithm  shown  in  Fig.  1  and  Fig.  2  gives  a  typical  structure  of  the  sawtooth 
noise  induced  by  the  receiver  and  provided  at  one  of  its  outputs  as  the  negative  sawtooth.  It  may  be  seen 
that  the  noise  is  mean-zero  with  the  upper  and  lower  bounds  of  ±50  ns  and  uniformly  distributed.  The 

A 

clock  first  state  estimate  x\ (n)  (TIE)  is  obtained  with  h^(i)  at  a  horizon  of  points.  The  observation 

A 

^2(ft)  for  the  second  state  x2(n)  (fractional  frequency  offset)  is  then  formed  by  increments  of  x\  ( n ) . 

A 

Accordingly,  x2{n)  is  achieved  with  h^i)  at  a  horizon  of  Nk~i  points.  Inherently,  the  first  accurate 

A 

value  of  x2{n)  appears  at  the  (NK  +  Nk~i) th  point  starting  from  ?z  =  0.  Finally,  the  last  state  estimate 

A 

Xk+ i  ( n )  is  calculated  with  h0(i)  at  a  horizon  of  No  points,  using  &+i(n)  that  is  formed  in  the  same  manner 

A 

as  ^2 (n).  The  first  correct  value  of  Xk+ i  (ft)  appears  thus  at  the  (NK  +  NK-i  +  . . .  +  Ao)th  point. 

The  unique  low  order  FIRs  for  the  algorithm  (Fig.  1)  are  derived  in  [4]  and  given  below.2 


K  (0  = 


N’ 


K(i) 


2  (2N  —  1)  —  6  / 
N  (TV  + 1) 


h2  (i) 


3  (3N2  -3N  +  2)  -18  (2AA-1)  /  +30/2 
7V(iv  +  l)C/V  +  2) 


(5) 

(6) 

(V) 


8  (2N3  - 3N2  +7N-3)-20  ( 6N 2  -6N  +  5)i  +120  ( 2N - 1)  i2  -  140  /3 

N  (N  +  i)  (N  +  2)  (A^  +  3) 


(8) 


having  the  inherent  common  properties: 


2  The  kernel  (6)  was  first  derived  in  [3]. 
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HO  =  fe»£r,  li*(0  =  1. 

i= 0 

Depending  on  the  model  degree  K ,  we  thus  have  several  particular  realizations  of  the  algorithm.  If  a 
noiseless  model  is  assumed  to  be  time-invariant,  K=  0,  the  only  nonzero  state  is  the  TIE  and  the  estimate 
of  the  TIE  is  provided  by  simple  averaging, 


1 


Nn-l 


^  0  i= 0 


The  linear  model,  K=  1,  is  processed  for  two  states  by 


x, 


i= 0 


(9) 


(10) 


Nn- 1 


A  1  JV0_A  A 

x2  in)  =  — —  (n  -  j )  - xx{n  -  j  —  1)]. 

For  the  quadratic  model,  K=  2,  the  3-state  unbiased  FIR  batch  algorithm  becomes 

a  n2-i 

xi  («)=  Xoc 

(=0 


W-l 


(2) = -  X  hi  0)[xi  (» -  y )  -  xi  o -  j  - 1)], 


(ii) 


(12) 


(13) 


A  ^  A 

X3(«)  =  —  —  ^[x2(n  -r)-jC2(«-r-l)]. 


rTV, 


0  r=0 


Finally,  if  we  assume  a  cubic  model,  K  =  3,  the  algorithm  is  formed  with 


A  ^3—1 

xi(n)=  EM'')  6  (»-*')> 

z=0 


(14) 


(15) 


A  J  X2-l 

x2  (»)  =  -  X  0  O')  [*i  0  “  i)  -  X\(n  -  j  - 1)], 


(16) 


A  J  Xl-1 

x3 0)  =  -  X h) O)  [*2 (n  —  p)  —  x2(n  —  p  —  1)],  (17) 

7”  ^=0 
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A  Y  N0  }  A  A 

x4(n)  =  — —  Y  [x3 (n - r) - x3 (n - r - 1)] ,  (18) 

tN0  % 

As  may  be  seen,  the  calculus  (9)— (1 8)  needs  setting  two  constants,  N  and  x,  for  each  of  the  estimates  to  be 
unbiased  and  near  optimal  in  the  sense  of  the  minimum  MSE. 


OPTIMUM  N  AND  x 

Experimental  evaluation  of  Nopt  and  Topt  for  each  of  the  unbiased  FIR  filters  in  (9)— (1 8)  was  provided  for 
the  crystal  clock  embedded  in  the  Stanford  Frequency  Counter  SR620  and  rubidium  clock  attached  to 
SR625.  Measurements  were  made  with  the  GPS  timing  sensor  SynPaQ  III  and  SR620  for  x  =  1  s  (GPS  - 
measurement).  Simultaneously,  to  obtain  reference  trends  for  each  of  the  estimates,  the  TIE  of  the  local 
clock  had  been  measured  for  the  atomic  clock  (we  exploit  the  rubidium  oscillator  attached  to  SR625  or 
cesium  clock  Cs  III  of  Symmetricom).  Initial  time  and  frequency  shifts  between  two  measurements  were 
then  eliminated  statistically  and  a  transition  to  x  >  1  s  was  provided  by  the  1PPS  measurements  thinning 
in  time.  The  difference  between  each  estimate  and  the  reference  trend  or  its  time  derivatives  is  evaluated 
in  the  sense  of  the  MSE  and  the  optimal  values  Aopt  and  iopt  are  found  for  the  minimum  MSEs. 

In  our  experiment,  for  signal  processing,  we  exploited  measurements  obtained  during  about  one  day. 
Therefore,  the  upper  bound  of  x  was  limited  by  104  s,  inducing  two  principal  limitations  for  Aopt  and  iopt: 


•  Finite  data  of  measurements  (about  1  day)  affecting  the  estimates  if  x  >  103  s. 

•  Long-term  phase  drifts  in  the  reference  clocks  affecting  estimates. 

Below  we  bring  results  of  experimental  evaluations  of  iopt  and  Aopt  for  several  states  of  a  local  crystal 
clock  via  the  GPS-based  and  rubidium-based  measurements. 

The  First  Clock  State  (TIE) 

Fig.  2  illustrates  an  evaluation  of  Aopt  vs.  x  for  the  first  clock  state  (TIE)  with  the  four  algorithms.  A 
common  conclusion  is  that  the  minimum  root-mean-square  error  (RMSE)  corresponds  to  x  =  1  s  by  all  of 
the  algorithms  and,  hence,  estimating  the  TIE  does  not  require  thinning  in  time  the  measured  data.  For 
TIE,  iopt  is  exactly  1  s.  One  may  also  observe  that  simple  averaging  provides  for  with  a  minimum  number 
Ao(opt)  of  the  points  among  all  other  filters.  However,  the  RMSE  in  simple  averaging  (5)  is  larger  than  in 
linear  and  quadratic  kernels,  (6)  and  (7),  respectively.  The  cubic  kernel  (8)  produces  even  larger  error 
than  in  (6)  and  (7).  Based  upon  this,  we  identify  the  TIE  model  of  the  investigated  clock  to  be  either 
linear  or  quadratic  on  the  horizon  of  Aopt  points. 

The  Second  Clock  State  (Fractional  Frequency  Offset) 

To  ascertain  iopt  and  7Vopt  for  the  second  clock  state  (fractional  frequency  offset),  one  first  needs  to 
estimate  optimally  the  TIE  for  the  optimal  N  and  x  taken  from  Fig.  2.  Thereafter,  by  the  algorithm  (Fig. 
1),  the  weighted  increments  of  these  estimates  serve  as  measurements  for  the  second  state.  Fig.  3  shows 
the  results  of  evaluating  the  Aopt  vs.  x  for  the  second  clock  state  and  we  indicate  that  the  minimum  RMSE 
corresponds  to  x  ranging  from  10  s  to  100  s.  We  notice  that  the  error  reduction  in  the  last  decade  of  x  in 
Fig.  3  is  caused  by  the  finite  database. 
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The  Third  Clock  State  (Linear  Fractional  Frequency  Drift  Rate) 

In  a  like  manner,  measurements  are  provided  for  the  third  clock  state  termed  the  linear  fractional 
frequency  drift  rate.  The  results  are  shown  in  Fig.  4.  Here,  a  minimum  RMSE  fits  the  region  of  x  from 
300  s  to  600  s.  For  this  minimum,  the  optimal  N  is  provided  by  Fig.  4  exactly,  even  though,  on  the  whole, 
the  Afyopt)  was  found  with  some  error.  We  notice  again  that  the  last  decade  in  the  x  scale  (from  103  s  to 
104  s  in  Fig.  4)  is  not  reliable  to  make  conclusions  for,  owing  to  the  finite  database. 

The  Fourth  Clock  State  (Quadratic  Fractional  Frequency  Drift  Rate) 

To  evaluate  Novt  and  Topt  for  the  fourth  clock  state  that  may  be  called  the  quadratic  fractional  frequency 
drift  rate,  the  only  algorithm  may  be  used,  that  is  (1 5)— (1 8).  Accordingly,  Fig.  5  gives  one  curve  for 
Afyopt)  and  the  other  one  for  the  RMSE  vs.  x.  The  principal  observation  is  that  the  minimum  RMSE  lies 
likely  in  the  last  decade  of  x  from  103  to  104  or  even  beyond  this  decade.  This  is  not  unexpected,  since  the 
frequency  changes  very  slowly  even  in  crystal  oscillators.  Therefore,  the  spectral  content  associated  with 
the  quadratic  drift  rate  inherently  places  very  closely  to  zero  and  much  time  is  necessary  to  evaluate  the 
Afyopt)  and  To(opt)-  Nevertheless,  one  may  conclude  that,  even  approximately,  the  quadratic  fractional 
frequency  drift  rate  may  be  estimated  by  the  unbiased  FIR  algorithm  with  Topt  about  3  hours  with  only  a 
couple  of  points  in  the  average. 


MEASUREMENTS  AND  ESTIMATIONS  OF  THE  TIE 

Provided  Topt  and  Nopt  for  every  clock  state  (Figs.  3-6),  we  compare  the  unbiased  FIR  estimates  to  the 
sawtooth  corrected  GPS-based  measurements. 

Almost  Linear  Trend  of  the  TIE 

In  the  first  experiment,  we  selected  the  part  of  the  process  in  which  the  TIE  of  a  crystal  clock  undergoes 
almost  linear  regular  changes.  Following  Fig.  3,  the  time  step  was  set  to  be  x  =  1  s  and  we  set  N=  2060 
and  N  =  2050  to  estimate  the  TIE  from  the  measurements  with  and  without  the  sawtooth  correction, 
respectively,  in  the  minimum  MSE  sense.  The  results  are  shown  in  Fig.  7.  The  first  point  to  notice  is  that 
the  RMSE  of  the  measurements  obtained  with  the  sawtooth  correction  is  4.67  ns  that  coincides  with  the 
facilities  of  the  correction.  The  linear  FIR  estimates  (6)  applied  to  measurements  without  and  with  the 
sawtooth  correction  were  found  to  be,  respectively,  1.6  ns  and  1.69  ns,  which  is  substantially  lower  than 
in  the  sawtooth-corrected  measurements.  An  additional  observation  shows  that  noise  in  the  estimates  is 
also  lowered,  having  no  excursions,  contrary  to  measurements.  Certainly,  for  the  near  linear  trend,  the 
unbiased  FIR  estimates  look  superior  to  the  sawtooth  correction. 

Almost  Quadratic  Trend  of  the  TIE 

For  the  second  experiment,  we  found  a  region  where  the  TIE  changes  almost  quadratically  (Fig.  8).  The 
optimum  numbers  were  found  here  to  be  N  =  920  and  N  =  790  for  measurements  with  and  without  the 
sawtooth  correction,  respectively.  We  notice  that  these  values  are  consistent  to  those  shown  in  Fig.  3  at  x 
=  1  s.  The  RMSE  error  of  the  sawtooth-corrected  measurements  was  calculated  to  be  7.5  ns.  Almost  the 
same  errors,  8.0  ns,  and  7.5  ns,  were  found  in  the  estimates  obtained  from  measurements  without  and  with 
the  sawtooth  correction,  respectively.  Because  the  nonlinear  case  is  more  regular,  we  conclude  that  the 
sawtooth  correction  and  the  estimation  produce  almost  the  same  RMSEs.  However,  like  in  Fig.  7,  the 
noise  in  the  estimates  is  still  substantially  lower  and  excursions  are  almost  fully  smoothed. 
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Complex  Function  of  the  TIE 

We  finally  examined  measurements  with  a  somewhat  sophisticated  behavior  of  the  TIE  (Fig.  9).  The 
optimum  numbers  were  found  to  be  N  =  1150  and  N  =  1120  (with  and  without  the  correction, 
respectively).  F  or  these  values,  the  RMSEs  of  the  estimates  were  calculated  to  be  7.83  ns,  and  8.0  ns, 
respectively,  and  the  RMSE  of  the  sawtooth-corrected  measurements  was  evaluated  by  8.87  ns.  Again, 
we  watch  for  the  comparable  results  and  notice  that  noise  and  excursions  are  still  much  more  suppressed 
in  the  estimates. 


MEASUREMENTS  AND  ESTIMATIONS  OF  THE  FRACTIONAL 
FREQUENCY  OFFSET 

By  the  algorithm  (Fig.  1),  measurements  and  estimates  of  the  second  clock  state  x2(n)  =  y  (n)  (fractional 
frequency  offset)  are  provided  with  the  discrete-time  derivatives  applied  to  the  TIE.  Accordingly, 

A  A  A 

measurements  of  the  fractional  frequency  offset  are  formed  by  \x\{n)  —  xx(n  —  1)]/t,  where  x\{n) 
should  be  treated  to  be  either  the  TIE  measurements  with  the  sawtooth  correction  or  the  estimates  of  the 
TIE  obtained  by  the  unbiased  FIR  filter.  By  the  results  given  in  Fig.  4b,  the  minimum  error  of  the  x2(n)  = 
y  ( n )  estimate  is  obtained  with  x  =  10  ...  100  s.  We  therefore  set  x  =  100  s,  since  it  corresponds  to  a  lower 
optimum  number  N  (Fig.  4a),  thus  a  smaller  computational  time.  The  results  of  measurements  and 
estimations  are  shown  in  Fig.  10.  Only  20  points  were  enough  to  provide  a  minimum  RMSE  in  the 
estimates.  We  therefore  set  N  =  20  to  a  simple  averaging  (5)  and  filter  the  time  derivatives  of  both 
estimates  of  the  TIE,  namely  with  and  without  the  sawtooth  correction. 

Fig.  10a  demonstrates  that  the  derivative  of  the  reference  measurement  fits,  on  the  whole,  the  derivative 
of  the  GPS -measurement  with  the  sawtooth  correction.  It  is  also  seen  that  the  GPS -measurement 
produces  much  larger  noise,  which  standard  deviation  reaches  55.6  x  10"12.  The  latter  value  may  be 
treated  as  a  resolution  in  frequency  of  the  GPS-based  sawtooth-less  measurements.  Contrarily,  the 
unbiased  FIR  estimates  (Fig.  10b)  produce  much  lower  noise  and  trace  very  close  to  Rb-measurements. 
Indeed,  the  standard  deviations  of  the  differences  between  the  reference  curve  and  two  estimates  (without 
and  with  the  sawtooth  correction)  were  evaluated  at  the  levels  of,  respectively,  15.0  x  10"12  and  13.4  x 
10"12.  Fig.  10c  demonstrates  how  well  the  estimates  fit  the  Rb-measurements.  For  the  comparative 
purposes,  in  Fig.  lOd,  we  show  the  estimates  obtained  by  the  linear  FIR  function.  For  this  filter,  the 
RMSEs  have  appeared  to  be  larger,  respectively,  16.7  x  10'12  and  14.6  x  10'12.  The  latter  means  that  the 
algorithm  is  correct  suggesting  the  uniform  kernel  (5)  for  the  second  clock  state. 


CONCLUDING  REMARKS 

In  this  paper,  we  presented  the  results  of  the  experimental  comparison  of  the  GPS-based  TIE 
measurements  (with  and  without  the  sawtooth  correction)  and  unbiased  FIR  estimation  for  a  local  crystal 
clock.  We  also  brought  the  relevant  plots  to  evaluate  the  optimal  time  step  Topt  and  optimum  number  of 
the  points  Aopt  in  the  average  for  the  uniform,  linear,  quadratic,  and  cubic  kernels  employed  by  the 
unbiased  FIR  filtering  algorithm  (Fig.  1).  The  conclusions  are  as  follows. 

The  optimal  values  Topt  and  Aopt  differ  for  different  clock  states.  For  the  crystal  clock  imbedded  in  the 
Stanford  Frequency  Counter  SR620,  the  first  state  (TIE)  must  be  filtered  with  x  =  1  s  and  the  relevant  Nopt 
taken  from  Fig.  3a.  The  second  state  (fractional  frequency  offset)  is  estimated  with  x  =  10  ...  100  s  and 
the  relevant  Aopt  provided  by  Fig.  4a.  The  third  state  (linear  fractional  frequency  drift  rate)  needs  setting  x 
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=  300  ...  600  s  with  Nopt  taken  from  Fig.  5a.  Finally,  the  fourth  state  (quadratic  fractional  frequency  drift 
rate)  may  be  evaluated  with  x  of  more  than  104  s  and  only  several  points  in  the  average. 

An  experimental  comparison  of  the  GPS-based  measurements  with  the  sawtooth  correction  and  the 
relevant  unbiased  FIR  estimates  of  the  TIE  model  of  a  local  crystal  clock  led  to  the  following 
conclusions: 

•  GPS-based  sawtooth-less  measurements  represent  basically  the  mean  value  of  its  origin  (without 
correction)  with  a  small  amount  of  noise.  Therefore,  the  unbiased  filter  applied  to  either  of  these 
measurements  produce  almost  the  same  output  (Figs.  7b,  8b,  9b,  and  10b). 

•  The  real-time  unbiased  FIR  filtering  algorithm  (Fig.  1)  efficiently  suppresses  noise  and  smoothes 
excursions  featured  to  sawtooth-less  measurements. 

•  The  unbiased  FIR  filter  produces  much  lower  errors  for  the  fractional  frequency  offset  and  higher 
degree  states  owing  to  the  lower  noise  intensity  in  the  estimates  of  the  TIE. 

Overall,  owing  to  small  noise  and  smoothed  excursions  in  the  estimates,  the  unbiased  FIR  filtering 
algorithm  (Fig.  1)  may  be  considered  to  be  more  efficient  than  the  sawtooth-less  measurements,  whenever 
accurate  evaluation  of  the  TIE  model  of  a  local  clock  is  needed. 
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Figure  1:  Structure  of  the  unbiased  FIR  filtering  algorithm  for  the  TIE  model  of  a  local  clock. 
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Figure  2:  Typical  structure  of  the  sawtooth  noise  induced  by  the  GPS  timing  sensor  SynPaQ  III. 


Figure  3:  Optimal  iVopt  vs.  r  for  the  first  clock  state  (TIE):  (a)  optimal  N  and  (b)  RMSE. 
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Figure  4;  Optimal  iVDpt  vs.  r  for  the  fractional  frequency  offset;  (a)  optimal  N  and  (b)  RMSJEL 


Figure  5;  Optimal  Nvp t  vs.  r  for  the  linear  fractional  frequency  drift  rate:  (a)  optimal  N  and  (b)  RMS'E, 


Figure  6:  Optimal  iV0pt  vs.  r  for  the  quadratic  fractional  frequency  drift  rate;  (a)  optimal  N  and  (b)  RMSE. 
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Figure  7:  Unbiased  FIR  filtering  vs.  the  sawtooth  correction:  (a)  GPS- based  nieasureiiienls  with  and  without 
the  sawtooth  correction  and  (b)  filtering  of  the  TIE  with  (light)  and  without  (bold)  the  sawtooth  correction 
using  the  linear  unbiased  FIR  filter.  White  curve  is  by  the  reference  Rb- measurements. 
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Figure3:  Unbiased  FIRfiltering  vs.  the  sawtooth  correction:  (a)  GPS-based  measurements  with  and  without 
the  sawtooth  correction  and  (b)  filtering  of  the  TEE  with  (light)  and  without  (bold)  the  sawtooth  correction 
using  the  linear  unbiased  FIR  filter.  White  curve  is  by  the  reference  Rb- measurements. 
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Figure  10'  Unbiased  FIR  filtering  of  the  second  dock  state  vs.  the  sawtooth  correction:  (a)  GPS- 
measurements  with  the  sawtooth  correction  and  Rb- measurements,  (b)  unbiased  FIR  filtering  of  the  mea¬ 
surements  with  (light)  and  without  (bold)  the  sawtooth  correction,  (c)  unbiased  FIR  filtering  with  a  uniform 
kernel  (5)  vs.  the  Rb- measurements,  and  (d)  unbiased  FIR  filtering  with  a  linear  kernel  (6)  vs.  the  Rb- 
measurements  (the  curve  labels  are  the  same  as  in  (c)). 
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